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Eliashberg theory is used to investigate the range of thermodynamic properties possible within 
a two-band model for s-wave superconductivity and to identify signatures of its two-band nature. 
We emphasize dimensionless BCS ratios (those for the energy gaps, the specific heat jump and the 
negative of its slope near Tc, the thermodynamic critical field Hc{0), and the normalized slopes of 
the critical field and the penetration depth near Tc), which are no longer universal even in weak 
coupling. We also give results for temperature-dependent quantities, such as the penetration depth 
and the energy gap. Results are presented both for microscopic parameters appropriate to MgB2 and 
for variations away from these. Strong coupling corrections are identified and found to be significant. 
Analytic formulas are provided which show the role played by the anisotropy in coupling in some 
special limits. Particular emphasis is placed on small interband coupling and on the opposite limit of 
no diagonal coupling. The effect of impurity scattering is considered, particularly for the interband 
case. 

PACS numbers: 74.20-z,74.70.Ad,74.25.Bt,74.25.Nf 



I. INTRODUCTION 

The properties of the superconducting state of conven- 
tional, single-band, electron-phonon superconductors dif- 
fer markedly from BCS predictions. 0| However, they are 
well-described within isotropic Eliashberg theory with a 
single electron-phonon spectral density a^F^ui) for the 
average interaction over the Fermi surface. This function 
is accurately known from inversion of tunneling data.|^ 
In many cases, the a'^F{uj) has also been calculated from 
first principle electronic band structure calculations ex- 
tended to include the electron-phonon interaction, some- 
times with the phonons taken directly from inelastic neu- 
tron scattering measurements. In many cases, such re- 
sults agree very well with the corresponding tunneling 
data. While it is to be noted that, in principle, the 
electron-phonon spectral density for the various electrons 
on the Fermi surface is anisotropic leading to energy gap 
anisotropy 0, 0, @1j this feature often does not play 
a prominent role because, in many instances, the elec- 
tronic mean free path is much smaller than the coherence 
length. In such circumstances, a Fermi surface average of 
the electron-phonon spectral density can be used. Nev- 
ertheless, corrections due to gap anisotropy have been 
I identified and studied in the past0, often, but not al- 
ways, in a separable anisotropic model. Q 

The history of two-band superconductivity^, 0, 

Sand of MgB, (with ~ 39 k|i3|) in particular^l 
El El Ei, 13 in m El is somewhat different. To 
our knowledge, as yet, there exists no inversion[2S| of 
tunneling data from which the electron-phonon interac- 
tion is determined. In fact, it has been noted |3| that 
this may well never be possible in MgB2 because of its 
two-band nature which requires a microscopic descrip- 
tion in terms of four separate electron-phonon spectral 



functions afjF{u!), where i = a,TT (or 1, 2), with the two- 
dimensional a band having the largest electron-phonon 
coupling. The three-dimensional tt band on its own would 
have a smaller value of Tc, the critical temperature, al- 
though it has a higher value of the electron density of 
states at the Fermi energy. 

In the absence of tunneling data giving reliable in- 
formation on the fundamental kernels entering the two- 
band Eliashberg equations, first principle band struc- 
ture calculations of afjF{uj) in MgB2 have been used 
to compute superconducting properties (for example. 
El El El E3)- To do this, it is also necessary to 
know the Coulomb pseudopotential repulsions /i*^ which 
are different for various indices but these have also 

been calculated. Good agreement with experiment is ob- 
tained in this way for the properties considered so far, 
more explicitly, the specific heat[l6ll25l|. the penetration 
depth [TtI l2^ . and the anisotropy in the two gaps, as 
well as their temperature dependence. For the penetra- 
tion depth, impurity scattering can be important, and 
in and out of plane orientation of the magnetic field are 
different Ell. 

In this paper, we use the band theory information on 
afjF{uj) and in MgB2 to calculate the critical tem- 
perature, the energy gap with its anisotropy and tem- 
perature dependence, and other thermodynamic prop- 
erties, as well as the penetration depth, giving partic- 
ular emphasis to strong coupling corrections. Further 
to our discussion of MgB2, we provide a full listing of 
calculated dimensionless BCS ratios, now modified by 
both the anisotropy and the strong coupling effects in 
MgB2, and make comparison with experiment. We also 
consider effects of variations in microscopic parameters 
away from those of MgB2 , as well as impurity scattering 
- intraband and interband. To this end, we reduce the 
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two-band Eliashberg equations, which fully account for 
retardation, in the two-square-well approximation (also 
called the A^^ model). This leads to simple renormalized 
BCS (RBCS) forms which, when compared to our full nu- 
merical Eliashberg results, allow us to identify the strong 
coupling corrections which we find to be significant even 
for MgB2. 

When considering variations in microscopic parameters 
away from those of MgB2 , we place particular emphasis 
on two limiting cases: the limit of small interband cou- 
pling and the opposite case, when the intraband coupling 
is zero and the superconductivity is due to the interband 
coupling alone, a case discussed in the early work of Shul 
et aL^9]. We also consider the special case when the in- 
traband coupling in the second band is repulsive. The 
limit of small interband coupling is particularly inter- 
esting because it allows us to understand how the off- 
diagonal terms lead to the integration of otherwise two 
completely independent and non-communicating super- 
conducting bands with separate transition temperatures 
Tci- In this regard, we find that a\2F{u)) and a\iF{Lo) be- 
have very differently with 21 the most effective variable 
at integrating the two systems and 12 the most effec- 
tive at changing the critical temperature. The presence 
of the offdiagonal interactions rapidly smear out the fea- 
tures of the second transition at Tc2, i.e, the one with the 
smaller of the two Td values. More specifically, surpris- 
ingly small values of the mass renormalization parameter 
A21, as compared with An and A22, have a large effect on 
the region of rc2. We also find that relatively modest val- 
ues of the interband impurity scattering rates lead to the 
significant integration of the two bands. Even when the 
bands are well-integrated, in the sense that little trace of 
a second sharp transition at Tc2 remains, there still exist 
important modifications of the usual one-band BCS re- 
sults because of the two distinct bands. As an example, 
the BCS dimensionless universal ratios now depend on 
the ratio of the electronic density of states at the Fermi 
energy of the two bands. Simple analytic expressions 
for these ratios are derived, which provide insight into 
the physics underlying two-band superconductivity and 
guidance as to how these results are to be interpreted. 

In section II, we give the two-band Eliashberg equa- 
tions and provide their reduction in the A^^ approxima- 
tion which is needed to identify strong coupling correc- 
tions to renormalized two-band BCS (RBCS). Section III 
deals with the dependence of on microscopic param- 
eters, ie, on the electron-phonon interaction as well as 
on impurities. Intra- and interband quantities are both 
of interest. We consider the modifications of the dimen- 
sionless BCS ratios in the A^^ model, as well as, the zero 
temperature value of the two gaps and their anisotropy. 
MgB2 is considered in section IV. The issue of strong cou- 
pling corrections in MgB2, and more generally in other 
related systems, is discussed. The limit of small inter- 
band electron-phonon coupling is considered in section 
V. We study, in particular, how the two otherwise sep- 
arate bands become integrated when this interaction is 



switched on. The effect of interband impurity scattering 
is also considered in the same context as it exhibits anal- 
ogous behaviour to the case of the offdiagonal electron- 
phonon coupling. In section VI, we deal briefly with the 
less realistic case of zero intraband electron-phonon cou- 
pling, where the superconductivity is due only to the in- 
terband piece. Conclusions are found in section VII. 

Finally, in light of the recent developments in other 
areas of superconductivity and correlated electrons, we 
wish to emphasize that our use of the term "gap 
anisotropy" here is in reference to the difference in 
the magnitudes of the two gaps, each of which are 
isotropic s-wave in this work, and hence does not refer 
to a momentum-dependent order parameter. Likewise, 
"strong coupling" refers to the traditional meaning of 
strong electron-phonon coupling and is not an allusion 
to strong interband coupling. 



II. THEORY 

The isotropic (within a band) Eliashberg equations 
generalized to two bands (i = 1,2) are written on the 



imaginary axis as 
Ai(iw„)Zi(iw„) = 7rTy^y^[Aij( 



and 



vrT 



Zi{iuJn) = IH '^'^^ij i^^rn " iw„) 



= , (2) 



where t+- = 1/{2ttt^-) and f^- = l/(27rr-^) are the ordi- 
nary and paramagnetic impurity scattering rates, respec- 
tively, and 



Qa^F.j (H) 



;dn. (3) 



Eq. ^ gives the gap Ai{iujn) and Eq. |(2Jl, the renor- 
malization Zi{iLUn) at the n'th Matsubara frequency iujm 
with ujn = (2rt — 1)ttT. Here, T is temperature and 
n = 0,±1,±2,---. The electron-phonon kernels are 
afjF{^l) as a function of phonon energy and the 
Coulomb repulsions are /z*^-, with a high energy cutoff 
LUc needed for convergence and usually taken to be about 



3 



six to ten times the maximum phonon frequency. For 
the specific case of MgB2, these may be found in [l^. 
The diagonal intraband elements of the electron-phonon 
interaction are largest, in the case of MgB2, while the 
offdiagonal elements describing interband scattering are 
smaller, but still substantial. 

In what is called the two-square- well approximation or 
A^^ model Jl, .27, 28], we use in Eq. P: 

\ij{iLOrn - ii^n) = , for both |w„|, \uJm\ < UJo 

= 0, otherwise, (4) 



where 

\ij {m ^ n) ^ Xi-j (0) = Ay 



dVL. (5) 



Neglecting the gap in the denominator on the right-hand 
side of Eq. (01 for Z, we further approximate (see Ref. ^ 
for details) 

Z,(iw„) = 1-l-^Ay. (6) 



This result may now be used in Eq. to obtain 

Ai(ztJ„) = Ai(T), \uJn\<UJo 

= 0, \uJn\ > OJo, (7) 



where 



A.(T) 



E E 

m,\ujm I <Wo j 



A,(T) 
4 + Af 



= [A.,-M*,], (8) 



where loq represents either the Debye frequency or 
some other characteristic energy scale representing the 
phonons in the system, at most the maximum phonon 
energy. Detailed justification of using a single cutoff is 
found in Ref. [ij. These results are used to derived 
various quantities within the A^^ model, which we will 
call renormalized BCS or RBCS. We also solve the full 
Eliashberg equations for typical strong coupling parame- 
ters and for the case of MgB2 , and in order to connect to 
the language most appropriate for this purpose, the mea- 
sure of the characteristic boson frequency, win, is defined 
to be: 



tJin = exp 



2 r°° 
— / Hcj) 

Ml Jo 



(9) 



This is reasonable for our case here as the wi„ calculated 
for the different afjF{Ljj) spectra of MgB2 are almost the 
same and other spectra used in this paper will have the 
same frequency distribution in each channel only scaled 
in magnitude. In general, this definition should be rea- 
sonably robust as, unless A22, A12, and A21, are large, the 
first channel An should dominate the strong coupling ef- 
fects. 



III. BCS RATIOS: THE A''" MODEL AND 
STRONG COUPLING 

A. Critical Temperature: Tc 

The critical temperature that results from the renor- 
malized BCS equation (|HJ| of the two-square- well approx- 
imation, takes the form 



A = In 



knTr 



or 



where 



(10) 



(11) 



^ _ All + A22 — \/ (All — A22)^ + 4A12A21 ^^^^ 
2(AiiA22 — A12A21) 



and 



T _ All - Mil T _ A12 - //12 

-^11 — :r~ — ^ — ' M2 — 



A22 — 



1 + All + A12 

A22 - M22 
1 + -^22 + A21 



1 + All + A12 

A21 = ~ ^^21 (^-^3^ 
1 + '^22 + '^21 



With no impurities and for one band (A12 = A21 = A22 = 
0) 



-l/Ai 



(14) 



Here we will be interested only in the ratio of Tc 
(Eq. inj) to (Eq. jTlJ)) and so the cutoff Wo can- 
cels, and the issue of the best choice for this quantity 
does not enter (see Allen and Dynes "2?!). Results for 
Tc/T°° based on Eqs. H11I14|I as a function of A21 for 
various values of A12 are shown in Fig. ^ where they 
are compared with results of complete numerical evalu- 
ation of the two-band Eliashberg equations (Q) and 
A Lorentzian model for the spectral densities afjF{uj) is 
used with zero Coulomb pseudopotential fi*^ for simplic- 
ity. Specifically, we use a truncated Lorentzian spectral 
density, which is defined in Ref. centered around 50 
meV with width 5 meV, truncated by 50 meV to either 
side of the central point. The win for this spectrum is 
44.6 meV. This spectral density is scaled in each of the 
four channels to give An = 1, A22 — 0.5, and the range 
of values of A12 and A21 as required for the figure. The 
curves, which are labelled in the figure caption, are for 
the renormalized BCS calculations and the correspond- 
ing Eliashberg calculations are presented as points. We 
note that for small values of A21 agreement between the 
A^^ results and full Eliashberg is excellent. The agree- 
ment is somewhat less good around A21 — 0.5 but still 
acceptable. An interesting point to note about this figure 
is that the effect on Tc of A21 and A12 are quite different. 
As A21 increases for fixed A12, Tc increases. On the other 
hand, for small but fixed A21, increasing A12 decreases Tc, 
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FIG. 1: Ratio of Tc to the pure, one-band T^" as a function of 
A21 for varying A12: 0.6 (long-dashed), 0.4 (short-dashed), 0.2 
(dotted), and 0.1 (sohd). Here, An = 1 and A22 = 0.5. Strong 
coupling Eliashberg calculations are given for comparison for 
the same parameters and are shown as the points with A12: 
0.6 (solid circles), 0.4 (solid triangles), 0.2 (solid squares), and 
0.1 (open circles). 



while the opposite behaviour is found to hold for values 
of A21 bigger than approximately 0.16. This behaviour 
is different from that expected in non-renormalized BCS 
theory where it is known that increasing the offdiagonal 
coupling from zero to some finite value always increases 
Tc whatever its sign. Expanding Eq. (|12|l under the as- 
sumption that the offdiagonal elements are small as com- 
pared with the diagonal ones (A12, A21 <C An — A22, A22) 
gives 



1 

1^7 



A12A21 
~1\ 



22 



1 



All - A 



22 



1 

A^T 



(15) 



In BCS theory, the Ay would not be renormalized as 
in Eq. H13|) . Since the term in curly brackets is posi- 
tive, A decreases with the product of A12A21 and hence 
Tc increases. But in our case, the multiplying term 
1/Aii = (1 + All + Ai2)/Aii contains A12 in leading or- 
der and this factor on its own increases A and there- 
fore decreases the critical temperature. These expecta- 
tions are confirmed in our full Eliashberg numerical work 
and are not captured in other BCS works (for example 
[30I l3]|). It is clear then, that in our theory, A12 and 
A21 do not enter the equation for in the same way 
because A12 provides a direct mass renormalization to 
the major interaction term An. If mass renormalization 
is ignored, as in BCS theory, this asymmetry no longer 
arises. The work by Mitrovic|32j on functional deriva- 
tives finds bTcl ^o?Fi\{J) to be positive and the one for 
12 to be negative, which conforms with our results. We 
note here that the disparity between A12 and A21, which 
will in turn affect the and other properties, is related 



to the different values of the density of states at the Fermi 
level Ni in each of the two bands, i.e. = N2/N1. 

Turning next to the effect of impurities on Tc, the 
change ATc = Tc — Tco for small impurity scattering can 
be written in the A^^ model as: 



AT, 



Tco All + A22 + 2A(Ai2A2i — A11A22) 



(16) 



where for ordinary impurities (C"*") and magnetic impu- 
rities (C~): 



= - — {(l-AA22)(p?2^1lTP^lAl2) 

+ (l-AAii)(p±A22T/0f2A2l) 

-I- AA2i(Ai2 T AiO/cifa 

+ ^Al2(A21 T A22)/5^i}, 



with 



P12 



cO 



1 + All + A 



\12 



P2I = 



cO 



1 + X 



22 



A 



21 



(17) 



(18) 



These equations have been derived for scattering across 
the bands; within the bands, paramagnetic impurities 
will affect Tc but ordinary, nonmagnetic ones will not. 
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FIG. 2: Ratio of Tc with impurity scattering to that with- 
out Tco as a function of tfJTco for varying A22: 0.5 (solid). 



All 



1, 



0.4 (short-dashed), and 0.3 (long-dashed). Here 
A12 = A21 = 0.02. For the lower three curves f^i = 
varies, and for the upper three curves (which are almost indis- 
tinguishable from each other) it is the reverse. In the middle 



and 



set of three curves, 



These calculations have been 



done with the full Eliashberg equations using a Lorentzian 
a'^F{uj) spectrum with Tco/i^'in ~ 0.11. The dotted lines are 
from the evaluation of Eq. (1161 for the A22 =0.5 case and are 
for = (upper dotted curve) and f J^ — (lower dotted 
curve). (Note that the middle set of curves are the only phys- 
ically realizable cases. The others serve to make the mathe- 
matical point that ™d ^21 affect Tc quite differently.) 
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Results are given in Fig.|5J Except for the dotted lines, 
all curves were obtained from numerical solutions of the 
linearized version of the Eliashberg equations Q and (0) 
using a Lorentzian model for af^F^uj). The curves come 
in sets of three for A22 = 0.5 (solid curve), 0.4 (dashed) 
and 0.3 (long-dashed). The other parameters are An = 1 
and A12 = A21 = 0.02 (small interband coupling). The 
lower set are for f^i = with varying while the up- 
per set have 1 <-> 2. The middle set have ~ ^21- 
Note that in showing the results when or i^j^ are var- 
ied separately, we are violating a requirement that they 
must be linked together by the density of states in the 
two bands. That is, as required for the Ay 's, likewise the 
impurity scattering rates must obey ^12/^21 = N2/N1. 
Our middle set of curves obey this constraint, but we 
have ignored it for the other curves in order to illustrate 
the general behaviour of each individual type of scatter- 
ing separately. As found for the A^-'s, the effect of t]^2 
and that of t^j^ on T^. are quite different. The quantity 
represents scattering from band 1 to band 2 and leads 
to pairbreaking much like paramagnetic impurities in the 
one-band case. We can see this analytically in the simple 
case of A12 = A21 = for which the two bands are de- 
coupled and the critical temperature is a property of the 
first band alone. In this case, Eas. H16ll8|l reduce to[3^ 



AT, 



T, 



cO 



-P12 



(19) 



for both normal or paramagnetic impurities in the linear 
approximation for the impurity scattering rate. The ini- 
tial linear decrease in with increasing p^2 is seen in the 
lower set of three curves of Fig. [3 As increased fur- 
ther, higher order corrections start to be important and 
the curves show saturation to a value which is larger, the 
greater the value of A22. Also note that formula H19|) 
shows that Tc is independent of p^i . This expectation is 
confirmed in the upper set of three curves of Fig. [3 where 
Tc has increased by no more than 3% for t+^/Tco = 1.5. 
This small increase is due to the small A12 = A21 used 
for the figure, while in Eq. H19|l . we have A12 = A21 = 0. 
The middle set of curves, which apply for — tti ^^"^ 
therefore satisfy the constraint imposed by having chosen 
A12 = A21 = 0.02, exhibits, by comparison to the other 
two cases, only a very small region which is linear in im- 
purity scattering and these curves are intermediate to the 
other two sets, as expected. They also saturate at higher 
values of Tc and we find that Tc decreases by only 20-30% 
for this case, similar to the observation by Mitrovic who 
was considering specifically the case of MgB2'3^. Fi- 
nally, we comment on the dotted curves which are based 
on Eqs. lfTB|) to lfTH|l valid in the A*^ model and first or- 
der in tfj. The lowest curve applies to the tj^ = case 

and the upper one to t^2 — 0. The slopes are in good 
agreement with the full Eliashberg results over a signifi- 
cant range of interband impurity scattering tf^ . For the 
middle set of curves the linear behaviour applies only 
comparatively to a rather small region. In all cases there 
still is some difference between A*^ results and Eliashberg 



because of strong coupling corrections. As previously 
stated, interband impurity scattering in two-band super- 
conductivity works like paramagnetic impurities in the 
ordinary one-band case. For this latter case, Schachinger, 
Daams, and Carbotte 35] have found for the specific case 
of Pb, the classic strong coupling material, that the A*^ 
model overestimates the initial slope of the drop in Tc 
value, with increasing impurity scattering. The physics 
is simple. For strong coupling, 2A/kBTc is larger than 
its BCS value i.e., the gap is bigger than expected on the 
basis of its Tc. This is because as T is increased, that part 
of the inelastic scattering which corresponds to the real 
(as opposed to virtual) processes, which are pairbreak- 
ing, increases and Tc is reduced below the value it would 
be without. As a result, the initial drop in Tc value with 
increasing impurity content is not as large in strong as in 
weak coupling because the system has a larger gap which 
is more robust against impurities. The same applies to 
interband scattering in a two-band superconductor. The 
initial slope of the drop is faster in the A^^ model than in 
Eliashberg, as most recently shown by Mitrovic 's^, who 
has commented on prior work by Golubov and Aiazin 33j , 
where only unrenormalized BCS results were given and 
the drop in Tc was even faster. Mitrovic also presents 
functional derivatives for ordinary impurities |34| and his 
findings compliment our calculations here. In addition, 
as low frequency phonons act like ordinary impurities, the 
previous work by Mitrovic on functional derivatives 32] 
for the electron-phonon spectral functions also confirms 
our impurity results by comparison with the behaviour 
of the low frequency part of the functional derivatives for 
12 versus 21. 

Finally, it has been of some interest amongst experi- 
mentalists, looking at novel superconductors, to know the 
outcome of having a repulsive interaction in the second 
band {i.e. A22 < 0). As will be seen in the next section, a 
second energy gap is still induced in this case due to the 
interband coupling, however, a signature of this repulsive 
band would exist in the case of impurity scattering, as 
strong interband scattering of sufficient strength could 
drive the Tc to zerO|3Qj. 



B. Energy Gaps and Gap Ratios 

We turn next to the consideration of the energy gaps. 
The transcendental equation for u = A2 / Ai at T = in 
the A^^ model is: 

— -^21 — — — — — — 

A12U h (A11A22 - A21A12) \nu = A22 - All, (20) 

u 

from which the gap ratio for the larger gap Ai may be 
found: 



In 



1.13Ai 

2kBTc 



A 



1 -|- Ai2Mlnu 



All + Ai2ii 



(21) 



The solution for the gap ratio 2Ai/fcB Tc can be corrected 
for strong coupling effects by multiplying by a factor ?7a 
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in the denominator of the logarithm of Eq. H21() with [36 



12.5( — 



In 



2Tr 



(22) 



As long as An is large and A22, A12, and A21 are small, 
one needs only to correct the first channel for strong cou- 
pling effects. Otherwise additional corrections for the 
other channels may exist but there would be no merit in 
such complexity of including these corrections over doing 
the full numerical calculations with the Eliashberg equa- 
tions. It is expected that in real systems, An is large rel- 
ative to the other parameters and hence dominates the 
strong coupling aspect of the result. However, when the 
offdiagonal couplings are significant, the strong coupling 
corrections of the first channel can affect the second. 
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FIG. 3: Gap ratios for the upper (2Ai/fc_B7c) and lower gap 
{2A2/kBTc) as a function of A12 for varying A22: 0.5 (solid), 
0.1 (dotted), -0.5 (short-dashed), and -1 (long-dashed). Here, 
All = 1, A21 = 0.3. These calculations are done using 
the RBCS formulas (I20I2H in the text, the solid dots show 
Eliashberg calculations for the same set of parameters with 
A22 =0.1 (for comparison with the dotted curve). Strong 
coupling corrections are significant and the rest of the curves 
in this figure would also be modified by strong coupling, much 
of this can be captured by the strong coupling correction for- 
mula given in the text. (Note that as A22 and A21 are finite, 
the points for A12 = are not physically realizable.) 

Our first set of results for the two energy gaps is given 
in Fig. 13 The lines are based on the simpler equations 
(|20(l and (|2I() . and the solid dots are for the results of 
full Eliashberg solutions on the imaginary axis and ana- 
lytically continued with Fade approximates ij to the real 
axis, where the gap is determined by Aq = A{uj = Ao)0]. 
For clarity in the figure, only one such set of results is 
shown for the case of A22 = 0.1. While magnitudes differ 
considerably between the renormalized BCS and strong 
coupling (comparing solid dots with the dotted curves). 



the general trends are the same. Specifically in Fig.O A12 
is varied with An = 1, A21 = 0.3, and A22 fixed to various 
values in turn. The upper curve applies to Ai and the 
lower curve of the same line type, to A2. While in all 
cases Ai increases with increasing A12, in one case (solid 
curve), the lower gap decreases slightly. More impor- 
tantly, the value of the upper gap ratio increases above 
its BCS ratio 3.53 and can reach 4.6 in renormalized BCS, 
a feature which comes from the two-band nature of the 
system. Comparing the dotted curves to the solid circles 
for Ai, we note that Eliashberg results are always above 
their A^^ counterpart, refiecting well-known strong cou- 
pling corrections to the gap. This applies as well to A2, 
the lower gap. We now comment specifically on the other 
curves. To increase the anisotropy between Ai and A2 
for the parameter set considered here, we need to de- 
crease the value of A22. Note, however, that even when 
we assume a repulsion in the second band, equal in size 
to the attraction An = 1 in the first band (long-dashed 
curve), a substantial gap is nevertheless induced in the 
second channel even for A12 = 0. It is the finite value of 
A21 which produces this gap. Recall that A21 describes 
the effect of band 1 on band 2 due to interband electron- 
phonon coupling. Turning on, as well, some A12 increases 
the second gap further but not by much. Finally, we men- 
tion that as A21 increases (not shown here), Ai decreases 
while A2 increases, ie. the ratio of A2/A1 goes up to- 
wards one and the anisotropy is reduced. 
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FIG. 4; Gap ratio 2A\ /ksTc as a function of A12 = A21, for 
All = 1.3 and A22 = 0.5. These curves provide a compari- 
son between the Eliashberg calculation (solid curve) and the 
RBCS calculation (dashed curve), along with the result from 
using the RBCS expression with the strong coupling correc- 
tion formula given in the text (dot-dashed curve). 

In Fig. 13 the ratio A12/A21 = N2/N1 is varying, while 
in Fig. ^ we keep A12 = A21 and illustrate more clearly 
the effect of strong coupling Eliashberg in comparison 
with the RBCS calculation, and also provide a compari- 
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son with the RBCS calculation corrected with the strong 
coupHng formula of Eq. One finds that the gap 

in Eliashberg is quite enhanced over the RBCS result, 
even exhibiting a different qualitative behaviour with the 
Eliashberg gap (solid curve) increasing with increasing 
offdiagonal A while the RBCS counterpart (dashed curve) 
is decreasing. However, when the strong coupling correc- 
tion formula is applied to the RBCS result, the resulting 
curve (dot-dashed) is now in reasonable agreement with 
the Eliashberg calculation and follows the evolution with 
increasing offdiagonal A very well. 

It is of interest to experimentalists '37] , looking at novel 
materials suspected of harbouring multiband supercon- 
ductivity, whether there may be a range of parameters 
that could produce a very large upper gap ratio with 
a large anisotropy in magnitude between the upper and 
lower gaps. It is possible that it could occur in a regime 
where A12/A21 » 1, as suggested by the trend in our 
Fig. 13 while in the opposite regime we will show that 
all results return to standard weak coupling BCS val- 
ues. As previously mentioned, this ratio of A12/A21 is 
equivalent to the ratio of density of states in the two 
bands, sometimes denoted as a in the literature, i.e. 
a = A12/A21 = N2/N1. We have gone to a = 20 within 
the renormalized BCS formalism and were not able to 
produce gap ratios bigger than about 5 or so, for the pa- 
rameters examined, and at the same time, the lower gap 
ratio was about 3. We conclude, therefore, that even with 
added strong coupling effects, very large gap ratios tend- 
ing towards 10 to 20 are difficult to obtain in conjunction 
with a large anisotropy in the two gaps. Repulsive poten- 
tials in the second band can give a large anisotropy, but 
they also lower the value of the upper gap ratio. Later 
in Section VI, we will return to this issue of trying to ob- 
tain large gap ratios and large gap anisotropy, when we 
examine another extreme limit first considered by Suhl 
et a/. 01. 

To conclude this subsection, we examine an ap- 
proximate formula for the gap ratio in two-band su- 
perconductivity, which has been given and used by 
experimentalists 38J , to determine its range of validity 
in the face of more exact calculations. The formula is an 
unrenormalized BCS formula and we have already seen 
that renormalization and strong coupling effects can be 
substantial. For A22,Ai2,A2i ^ An, we can derive the 
primary (or large) gap ratio as: 



only constraint being that u ^ 1. 



2Ai 



3.53 



3.53 



^^12 2 1 

- — u In-i 
A21 



7V2M2V, f^2 



(23) 



which is the same equation as given in lavarone et al. j38l | , 
where their use of the indices 1 and 2 are reversed with 
respect to ours. In our formula H23|l given here, the u 
and A's are coupled through Eq. ((SOJ, but in the case of 
Ref. (3^ the ratio of the density of states and the ratio of 
the gaps are treated as independent parameters with the 
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FIG. 5: Gap ratios for tlie upper {2Ai /kBTc) and lower 
gap {2/S.2/kBTc) as a function of A12 for An = 1.0, A22 = 
0.5, A21 = 0.2. The solid curve is the exact BCS result, 
whereas, the dashed curve illustrates the approximate formula 
of lavarone et al. 1331 . 



In Fig. [SI we compare this approximate BCS formula 
with that of our exact renormalized BCS formula for typ- 
ical \ij values used in the literature. The /i*^ are set to 
zero as there is no such feature in the lavarone et al. 
formula and the /Lt*'s in that case would simply serve to 
change the effective value of A's. We find that the ap- 
proximate formula (dashed curve of Fig[SJ) compares well 
with the renormalized BCS result in the limit of small 
Ai2,2i.22, as required by the constraint of the approxima- 
tion, and breaks down for A12 > 0.5, where the approxi- 
mate formula tends to overestimate quite significantly the 
value of the two gaps. Strong coupling effects would pro- 
duce very significant deviations in addition. Not shown 
is the case where Ai2,2i,22 were all taken to be very small 
and then in that expected, there was excellent 

agreement between the exact renormalized BCS calcula- 
tion and the approximate form. The fact that lavarone 
et aLjs^ obtained excellent estimates of the two energy 
gaps for MgB2 is maybe fortuitous in some sense, because 
it will be seen in the next section, where we discuss MgB2 
in detail, that the renormalized BCS formula underesti- 
mates the correct gap values of MgB2 and strong coupling 
corrections of about 7-10% are needed to obtain good 
agreement between the data and full Eliashberg calcula- 
tions. We conclude that their simple formula is helpful, 
but that it should be used with caution when consider- 
ing systems where the parameters are no longer small as 
then this formula will fail. 
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C. Specific Heat Jump 

The specific heat is calculated from the free energy. 
The difference in free energy AF = Fs — F^ between 
the superconducting state and the normal state is given 
byQ: 



+ 00 



n— — 00 t 



(24) 



where "S" and "N" refer to the superconducting and 
normal state, respectively, and i indexes the number of 
bands. From this, the difference in the specific heat is 
obtained: 



AC = -T- 



^£AF 



(25) 



and the negative of the slope of the difference in specific 
heat near Tc is given as 



dAC(T) 



dT 



1 



(26) 



where 7 is the Sommerfeld constant for the two-band 
case. 

In the A®^ model, the gap near Tc, for t — T/Tc, can 
be written as 



A?(t) = ^^(l-t), 
7C(3) XI 



(27) 
(28) 



7C(3) X2 
where C(3) ~ 1.202. Here, 

(1 ^ AA22)Aii + Ai2A2iA[l + A^'Xj^jl ~ A\^2Y^\ 

(1 - ^A22)Aii + Ai2A2i[2A + ^2X22/(1 - AA22)] 

(29) 



Xi = 



and 



X2 



(1 - ^An)A22 + A2iAi2A[l + A^\\^(\ - AAii)-^] 
(1 - AAn)A22 + 2AA2iAi2[2A + A^\xl{\ - AAn)] ' 

(30) 

and the strong coupling correction is introduced 
through js^: 



(31) 



The specific heat jump at Tc is 
AC 



= 1.43 



(1 + All + Ai2)^ + a{l + A22 + A2i)^ 



(1 + All + A12) + a{l + A22 + A21) 



(32) 



We find with this expression that anisotropy (ie., An 7^ 
A22) reduces the jump ratio but increasing A12 or A21 
increases the ratio, and the maximum obtainable is 1.43. 
Other work along the same line is given in Refs. 30, 
i31j where they do not consider full renormalized BCS or 
strong coupling theories, as we have done here. 

When Ai2,A2i ^ 0, 1/xi 1 + 0(\l^) and l/x2 ~ 
0(^12) • This is assuming An — A22 and A22 remain sig- 
nificant as compared with the value of the offdiagonal 
elements. In this case. 



AC 



= 1.43 



(1 All A12) 



(1 + All + A12) + a{l + A22 + A21) 



(33) 



The physics of this formula is that, in this limit, the 
specific heat jump at Tc itself is determined only by the 
superconductivity of the dominant band, but it is nor- 
malized with the normal state specific heat 7 belonging 
to the sum of both bands. This has the effect of making 
AC{Tc)/jTc always less than the BCS value by a factor 
of l/(l-ha*), where a* = a(l-hA22 + A2i)/(l-hAii-l-Ai2). 

For MgB2, we expect a* ^ 1 which means that in this 
case the normalized jump is reduced to about half its BCS 
value. If we had included in (|33|l the strong coupling 
correction rjc, this would have the effect of increasing 
the factor 1.43 to a larger value characteristic of strong 
coupling but the additional anisotropy parameters would 
still work to reduce the jump. Thus, in a two-band su- 
perconductor, the jump will be smaller than for one band 
with the same strong coupling index 



D. Thermodynamic Critical Magnetic Field 

The thermodynamic critical magnetic field is calcu- 
lated from the free energy difference: 



Hc{T) = V-SvrAF. 



(34) 



As the temperature dependence of this quantity, normal- 
ized to its zero temperature value, follows very closely a 
nearly quadratic behaviour, the deviation function D(t) 
is often plotted: 



where t ^ T/T^ 
At T = 



(35) 



(36) 



i/2(o) ^ 47riVi*A2(l + a*u2), 

where a* = N2/NI and 

iV* -7V,(0)(1 + A,, + A,,). (37) 

The zero temperature critical magnetic field is modified 
through the second term in H36|) which increases with in- 
creasing a* and with the square of the anisotropy ratio 
u, which in this case is just the ratio of the independent 
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gap values for the two separate bands. Further, the di- 
mensionless ratio is 



_ 7r(fcsT,)2[i + a*] 



772(0) 6A?[l + a*ii2] 
For almost decoupled bands, Eq. H38|) becomes 



(38) 



0.168, 



1 



1 + a*u'^ 



(39) 



where the second factor on the right-hand side modifies 
the usual single-band BCS value of 0.168 for the presence 
of the second band. Again, both a* and u enter the 
correction. If there is no anisotropy, u = 1, and therefore 
the bands must be the same, we recover the one-band 
limiting value. For large anisotropy where m — > 0, and if 
a* is of order one, the ratio in Eq. H39|l is of order twice 
its one-band value because the second band contributes 
very little to the zero temperature condensation energy, 
but is still as equally important as the first band in its 
contribution to ^Tc, the normal state specific heat. Near 



/ 327r 
703) 



(^fcBT,)(l-i) 



1/2 



(40) 



which then gives the dimensionless ratio 



\mTc)\T. 

2Ai 1 



ksTc TT V 32 



7C(3) 



/ l + a*u2 

-2 , I ~' 

Xl +a*X2 



(41) 



Strong coupling factors could be introduced in H36|l . 
(|38|l . and (|40|) . They are not given explicitly here as 
they are less important than for the specific heat jump 
and the slope of the penetration depth at Tc (see Ta- 
ble P). The limit of nearly decoupled bands (A12, A21 <C 
All — A22, A22) gives for this quantity: 



hc{Q) = 0.576VT 



-I- a*u^. 



(42) 



The square root, which accounts for two-band effects con- 
tains a correction proportional to a*v?'. It can be un- 
derstood as follows. The slope at Tc found from formula 
(|40|l depends only on band 1 but i7c(0) involves both and 
hence this correction comes solely from Hc(0) as seen in 
Eq. H36|l . If the anisotropy between the two bands is large 
M — > 0, there is no correction factor in (|42|l because the 
second band is eliminated from Hc{0). If, on the other 
hand, u is near 1, the two bands have nearly equal gap 
value but still it is only band 1 which contributes to the 
slope at Tc and the dimensionless ratio (|42|) can now be 
larger than its BCS value. 



E. Penetration Depth 

The London penetration depth Xl{T) is evaluated 
fromQ: 



1 



Ai(r) 



Af(iw„) 



(43) 

where in three dimensions 



1 



N,v] 



(44) 



and vpi is the Fermi velocity in the band labelled by 
the index i. This last equation would be multiplied by a 
factor of 3/2 in two dimensions. 

For the penetration depth Al(T) at T = 0, 



Ai(0) vl.mioii^ + An + Ai2)^A2„2(l + A22 + A2i)^ 

(45) 

and near T^, 



Ai(t) 



2(1 -t) 



1 



<(7^c)ALiXi(l + Aii + Ai2) 
1 



where 



AL2X2(1 + A22+A2l) 



(46) 

(47) 
(48) 



Hence, defining yL{T) — 1/X\[T), we write the dimen- 
sionless BCS penetration depth ratio y as 



=1(1 



-a/3) 



1 

Xl 



a(3 
Xl 



(49) 



where (i = v\^{\ + An + \\i)lv\^{\ + A22 + A21). /3 is 
expected to be of order 1 unless there is a great disparity 
in the two Fermi velocities. For MgB2, we use the values 
of vpx = 4.40 X 10^ m/s and vp2 = 5.35 x 10^ m/s re- 
ported in Ref. ^^'^ for our other model calculations, 
we take them to be equivalent, for simplicity. For the 
nearly decoupled case 



1 



\v'L<rc)\Tc 2 



(l + a/3). 



(50) 



For a and /3 equal to one, we see that the normalized 
slope of the penetration depth is twice its one-band BCS 
value of 1/2. Should a, /3, or both be much larger than 
1, then the slope can be even larger, which reflects the 
fact that only the dominant band determines the slope 
y'j^ but both bands contribute to ?/l(0). Information on 
the vpi and the iVi(O) is contained in the slope. 
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FIG. 6: Upper frame: Electronic specific heat for MgB2 in 
the superconducting state normalized to the normal state 
as a function of T/Tc. The points are the experimental re- 
sults of Wang et al. |25|I and the solid curve is the result for 
the Eliashberg calculation using the parameters given in the 
literature |3^. The dashed curve illustrates the case where 
the Ai2 and A21 parameters, used for the solid curve, have 
been halved. The jump due to the lower gap begins to ap- 
pear in this case. Middle frame: [A(0)/A(T)]^ versus T/Tc. 
Curves are those resulting from the same set of parameters 
as discussed for the upper frame, with the vpi taken from 
Ref. [1^. The data, shown for comparison, have been taken 
from Ref. [l^ . No impurity scattering has been used to ob- 
tain a better fit. Lower frame: The deviation function D{t) 
for the thermodynamic critical field. Line labels are as above 
and the data (open and solid circles) are formed from the 
Hc{T) data given by Wang et al. |25| and Bouquet et al. |40( . 
respectively. 



IV. MgBa: INTEGRATED BANDS AND 
STRONG COUPLING 

We now continue on beyond renormalized BCS for- 
mulas to evaluate quantities based on the full two-band 
Eliashberg formalism and we begin with the specific 
case of MgB2 and strong coupling effects. Eqs. and 
(0) were solved for electron-phonon spectral densities 
a-fjF{uj), read from graphs in Ref. '35|, which were orig- 
inally presented in Ref. 19]. The Coulomb repulsion pa- 
rameters iJ,*j and Xij, taken from were: Ao-o- = 1.017, 
A^^ = 0.448, A,,^ = 0.213, = 0.155, = 0.210, 
fil^ = 0.172, = 0.095, and /.t;^ = 0.069, with 
LUc — 750 meV. From these parameters, Tc was found 
to be 39.5 K. As discussed in our theory introduction, we 



used Win = 66.4 meV, calculated from the aliF{uj) spec- 
trum, to form our strong coupling index Tc/uj\n- The 
other three channels had win — 62 meV, which is not so 
different, although as argued previously, the main strong 
coupling effects will come from the 11 channel, and hence 
the choice of 66.4 meV for this parameter. From the solu- 
tion of the Eliashberg equations, we can evaluate Eq. H24|) 
for the free energy difference between the superconduct- 
ing and normal state, and evaluate the superfluid den- 
sity or the inverse square of the penetration depth from 
Eq. H43|) . In Fig.|H| which has three frames: the top is the 
specific heat, middle, the penetration depth, and bottom, 
the critical magnetic field deviation function of formula 
(|35ll , we compare Eliashberg results (solid curve) with ex- 
perimental results (solid and open circles, triangles, and 
squares). 

In all cases, the agreement with experiment is very 
good and certainly as good as is obtained in conventional 
one-band cases 0. In each case, we also present a sec- 
ond set of theoretical results (dashed curve) for which 
all microscopic parameters remain those of MgB2 except 
that we have half the value of the offdiagonal spectral 
functions a\2F{Lo) and a\iF{uj), which changes the Tc 
only by about one degree. It is clear that doing this 
reduces greatly the quality of the fit one obtains with 
the experimental data. This can be taken as evidence 
that the electronic structure, first-principle calculations 
of electron-phonon spectral functions are accurate. It 
also shows that variation of parameters by a factor of 
two or so away from the computed ones can lead to sig- 
nificant changes in superconducting properties and, in 
this instance, features of the second transition, due to 
the lower gap, begin to appear. The specific heat curve 
was computed before in Refs. |0, and the penetra- 
tion depth in Refs. [T^ [26|. In these cases, our calcu- 
lations (solid curves) confirm previous ones and demon- 
strate that our calculational procedure is working cor- 
rectly. For the penetration depth we did not introduce 
impurity scattering. Impurities can affect the penetra- 
tion depth and were included in Ref. |0 . The three sets 
of penetration depth data are for clean (solid circle s ISOj 
and triangles li^) a nd dirty samples (solid squares [Slf) 
as discussed in T7]. To our knowledge, the deviation 
function has not been computed and compared with ex- 
periment before and it is presented for the first time here. 
The data is from Refs. [23 (open circles) and 40] (solid 
circles) and again agreement with calculation, with no 
free parameters, is very good. The minimum in the de- 
viation function for the Eliashberg calculation occurs at 
T/Tc = 0.6 and has a value of -0.054. In the experimental 
data, the minima occur at about T/Tc = 0.6 and 0.65, 
with values of about -0.05 and -0.045, respectively. For 
reference, the one-band BCS value is -0.037 and strong 
coupling makes this value even smaller and can even push 
it to a positive value, hence anisotropy is compensating 
for the strong coupling effects and is making this value 
larger and more negative. 0] 

In Fig. [71 we present the temperature dependence of 
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FIG. 7: Gap ratios for the upper (2Ai/A;bTc) and lower gap 
{21S.2 /knTc) as a function of T jT^. Shown as the solid curves 
are the predictions for the gap ratios given by our full Eliash- 
berg calculations for MgB2, the dashed curves are the Eliash- 
berg calculations for the case of reducing the offdiagonal A's 
by half and the dotted curves show the classic BCS temper- 
ature dependences to illustrate the deviation of the temper- 
ature dependence of the Eliashberg two-band calculation for 
MgB2. The open circles are the data from lavarone et al. |38l . 
where we have used a Tc = 38.3K to obtain their quoted 
upper gap ratio value of 4.3. The solid dots are the data of 
Gonnelh et al.^. 



the two gap ratios for MgB2. Once again the sohd curve 
is the full Eliashberg calculation using the parameters 
given for MgB2 with no adjustments. The ratio A1/A2 
increases from 2.7 at T = to about 3.5 at Tc- The 
temperature-dependent behaviour shown here was also 
found by Choi et a/.0], Brinkman et al.^^, and Gol- 
ubov et aL|l9l|. A comparison with some of the more 
recent experiments is given by the open and closed cir- 
cles, with the data taken from lavarone et a/.|38l| and 
Gonnelli et al. l4l| . respectively. Similar data is found in 
other references [43. \4M . l45l |. In the case of the data by 
lavarone et al., the statement of Tc was ambiguous and 
so we used their quoted value of the upper gap ratio of 
4.3 along with their quoted value of the upper gap being 
7.1 meV to determine a Tc = 38. 3i^ used for the scaling 
of the data for the plot presented here. The Gonnelli et 
al. data is presented based on the Tc of 38. 2K given in 
their paper. There is a very reasonable agreement of the 
data with the calculation, once again, along with Fig. El 
this shows a consistency of a number of sets of data from 
several different types of experiments with the one set of 
parameters fixed from band structure for MgB2. Thus 
overall, the agreement between theory and experiment is 
excellent and validates the two-band nature of supercon- 
ductivity in this material. The dotted curves in Fig. [7| 
are presented to show that the two-band calculations do 
show deviation from a classic BCS temperature depen- 



dence (which was used in the original presentations of the 
datajs^^). In particular, Gonnelh et al. argued that 
the deviation of their lower gap data at temperatures 
above 25K (or T/Tc = 0.65, here) from the BCS temper- 
ature dependence is an additional signature of the two- 
band nature of the material. However, we find no such 
dramatic suppression in the two-band calculations at this 
temperature and only with the dashed curve, where we 
have taken the offdiagonal electron-phonon coupling to 
be half of the usual value for MgB2 do we find an in- 
fiection point around 0.35. We were not able to induce a 
suppression of the lower gap in the vicinity of Tc by vary- 
ing the MgB2 parameters slightly about their accepted 
values. However, such behaviour can be found in other 
regimes of the parameter space not relevant to MgB2 and 
this feature and the issue raised by Gonnelli et al. will be 
discussed further in the next section. To end, note that 
an inflection point is also seen in the penetration depth 
at about T/Tc ~ 0.35, as described first by Golubov et 
al.\n^ and also found here (solid curve of middle frame 
of Fig.EJl. 

More results from our calculations as well as compar- 
ison with data are presented in Tabled In the first col- 
umn, we include, for comparison, the one-band BCS val- 
ues for the various dimensionless ratios. The strong cou- 
pling index is first, followed by the major gap to critical 
temperature ratio, the minor gap ratio, the anisotropy 
A2/A1, the normalized specific heat jump and the neg- 
ative of its slope at Tc, jT^/H^{0), and the inverse of 
the normalized slope at Tc for the critical magnetic field 
and for the penetration depth. Included in the second 
column, also for comparison, are the same indices for 
Pb, the prototype, single-band, strong coupler. We re- 
mind the reader that, in many conventional superconduc- 
tors, strong coupling corrections are large and that the 
data cannot be understood without introducing them, 
and these are to be differentiated from those corrections 
due to anisotropy. The third column gives the results of 
our two-band calculations for MgB2. This is followed by 
a column giving experimental values. It is clear that the 
agreement between theory and experiment is good. Note 
that we have not attempted to make a complete survey 
of all experiments, but have tried to present as many 
as reasonable, with no judgement about the quality of 
the data or samples, which might have improved over 
time. In addition, for the quantities related to slopes, 
i.e., g, hc{0), and y, we have tried to estimate these our- 
selves from the graphs in papers and so this should be 
viewed as rough estimates as the values might change 
with a more rigorous analysis of the original data. Also 
shown are the results when our renormalized BCS formu- 
las of the previous section are implemented using MgB2 
parameters 52J , which allows us to define a measure of 
strong coupling corrections, entered in column 6 as per- 
centages. It is seen that MgB2 is an intermediate cou- 
pling case. The next column shows the results when the 
analytical expressions for strong coupling corrections to 
renormalized BCS, given in the text, are applied. This 
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TABLE L Universal dimensionless BCS ratios and their modification for strong coupling (SC) and two-band superconductivity. 
RBCS stands for Renormalized BCS formula given in text. The percentage difference between the full Eliashberg (Eliash.) 
calculation and RBCS, used to measure the amount of strong coupling correction, is given as % SC and defined as | (Eliash. — 
RBCS) /Eliash. I . 



Ratio 


BCS 


Pb 


MgBa 


MgBa 


MgBa 


MgBa 


MgBa 


Lor 


Lor 


Lor 


Lor 




one-band 


one-band 


Eliash. 


Expt. 


RBCS 


% SC 


RBCS-HSC 


Eliash. 


RBCS 


% SC 


RBCS-hSC 


Tc/uJln 


0.0 


0.128 


0.051 


0.076" 


0.0 




0.052 


0.15 


0.0 




0.15 


2Ai/kBTc 


3.53 


4.49 


4.17 


3.6-4.6' 


3.86 


7.4% 


4.15 


4.97 


3.84 


23% 


5.14 


2A2/kBn 


3.53 


4.49 


1.55 


1.0-1.9' 


1.40 


9.7% 




2.66 


2.27 


15% 




A2/A1 


1.00 


1.00 


0.37 


0.30-0.42' 


0.36 


2.7% 




0.535 


0.593 


11% 






1.43 


2.79 


1.04 


0.82-1.32" 


0.817 


21% 


1.02 


2.08 


1.07 


49% 


1.97 


9 


-3.77 


-12.68 


-3.28 


-(2.37-4.31)'' 








-8.32 
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improves the agreement vifith the full Eliashberg results 
as compared to RBCS. Some discrepancies remain due 
in part to additional modifications introduced by the 
coupling of a strong coupling band with a weak cou- 
pling one through the offdiagonal Ay's. The next four 
columns were obtained for our Lorentzian spectral den- 
sity model with An = 1.3, A22 = 0.5, A12 = A21 — 0.2, 
and = 0. This was devised to have a strong cou- 
pling index Tc/luij^ ~ 0.15 which is slightly larger than Pb 
and well within the range of realistic values for electron- 
phonon superconductors. It is clear that strong coupling 
corrections are now even more significant and cannot be 
ignored in a complete theory. 

More information on strong coupling effects as well as 
on two-band anisotropy is given in Fig.|Hl where we show 
the same BCS ratios as considered in Tabled In all eight 
frames, we have used our model Lorentzian afjF{u!) spec- 
tra. The solid curves are results of full Eliashberg calcu- 
lations as a function of A12 = A21, with An fixed at 1.3 
and A22 at 0.5. The dashed curves are for comparison 
and are based on our A^* formulas, i.e, give renormalized 
BCS results without use of the strong couphng correction 
formulas. They, of course, can differ very significantly 
from one-band universal BCS values because of the two- 
band anisotropy. We see that these effects can be large 
and that on comparison between the solid and dashed 
curves, the strong coupling effects can also be signifi- 
cant. As A12 = A21 is increased from zero, with An and 
A22 remaining fixed, the Tc increases and this leads to 
the increase in Tc/win from about 0.15 at A12 = A21 — 
to over 0.2 at A12 = A21 = 1. For all the indices consid- 
ered here, we note that their values at Tc/win = 0.2 are 
close to the values that they would have in a one-band 



caseQ, and the remaining anisotropy in the Ay's play 
only a minor role. (Of course, this is a qualitative state- 
ment since it is well known that the shape of a^F(cj) for 
fixed T^/win can also affect somewhat the value of BCS 
ratios0|.) This is expected since in this case the fluctu- 
ation off the average of any Ay is becoming smaller. For 
RBCS all ratios have returned to the one-band case at 
A12 = A21 — 1 except for y which remains 6% larger. We 
now comment on select indices separately. The normal- 
ized specific heat jump at Tc in the A^^ model is given 
by formula H32() with r/c = 1- For A12 — X21 small, 
Xi^ ^ 1 + Om^) and ^2"^ - + OiXj^). These con- 
ditions mean that AC/jTc rises slightly as A12 = A21 
increases, and eventually reaches 1.43. By contrast, the 
solid curve includes, in addition, strong coupling effects 
which increase the value of the jump ratio rather rapidly. 
For 2Ai^2/fcB?c, the lower gaps have the same value for 
A12 — A21 = as it is determined only by A22- This is 
not so for the upper gaps. The dashed curve takes on 
its BCS value of 3.53, but the solid curve (an Eliashberg 
calculation) has strong coupling effects as described in 
Fig. El (This means that A2/A1 is smaller for the solid 
curve as compared to the dashed one in the lower left- 
hand frame.) As A12 = A21 increases, the long-dashed 
and lower short-dashed curves begin to deviate because 
the former starts to acquire strong coupling corrections of 
its own through the offdiagonal A's. While the solid curve 
also increases, the anisotropy between 1 and 2 decreases. 
The short-dashed curves show different behaviour. The 
ratio 2Ai/fcsT'c starts at 3.53, rises slightly towards 4 
before tending towards 3.53 again. Now, the anisotropy 
between A2 and Ai decreases mainly because A2 itself 
rises towards 3.53. The behaviour of 'yT^/H^{0) (dashed 



13 



u 
<1 



0.3 



0.2 



0.1 
3 



1 


- 0.8 
< 0.4 



0.0 



-8 

Cifl -10 

-12 

-14 
0.3 




0.2 

0.1 
^ 0.7 

o 

^ 0.6 
0.5 
0.4 
1.2 
^ 0.8 
0.4 



0.0 



0.5 

A 



1.0 



0.0 




0.0 



12 



0.5 

A 



1.0 



12 



FIG. 8: Various BCS ratios as discussed in the text, shown as 
a function of A12, where A21 = A12 (i.e. a = 1), An = 1.3, and 
A22 =0.5. The solid curves are those for the full Eliashberg 
calculation for a Lorentzian model of o?Fij {uj) spectra and the 
short-dashed curves are for the renormalized BCS formulas 
developed from the A*" model and given in the text. For the 
frame with the gap ratios, the upper gap is given by the solid 
curve and the lower gap is given by the long-dashed curve, 
the upper and lower short-dashed curves are for the upper 
and lower gaps, respectively, in RBCS. The first frame gives 
the effective Tc/uii-n for the Eliashberg spectrum based on the 
definition given in the text. 
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FIG. 9: Upper frame: Specific heat in the superconduct- 
ing state normalized to the normal state, CsiT)/yT, versus 
T/Tco, where Tco is the Tc for only the An channel, with all 
others zero. Shown are curves for various offdiagonal A's with 
All = 1 and A22 = 0.5. Three curves are for A12 = A21 equal 
to: 0.0001 (solid), 0.01 (short-dashed) and 0.1 (long-dashed). 
Also shown are: A12 = 0.1 and A21 ~ 0.01 (i.e., a = 10) 
(dot-dashed) and A12 = 0.01 and A21 = 0.1 (i.e., a = 0.1) 
(dotted). Middle frame: The superfluid density [A(0)/A(T)]^ 
versus T/T^o for the same parameters. Lower frame: The de- 
viation function D{t) plotted versus T/Tco- The dot-dashed 
curve has been divided by 10 from its original value in order 
to display it on the same scale as the other curves. 



curve) can be understood from Eq. (|38|l . While Ai/Tc, 
as we have seen, does change somewhat with A12 = A21, 
a more important change is the factor in the denom- 
inator of (|38|l which rapidly decreases this ratio towards 
its BCS value of 0.168 as u increases towards 1. The 
behaviour of hc{Q) given by Eq. (|41|l is more complex. 
The numerator in the square root goes towards 1 + a* , 
as — > 1, more rapidly than the denominator which 
involves the x's- Here, the numerator and denomina- 
tor compete and consequently hc{0) first increases before 
showing a slow decrease to its BCS value. Finally, y in 
formula H49|) decreases with increasing offdiagonal A be- 
cause of the square bracket in the denominator. It is clear 
from these comparisons between Eliashberg and RBCS 
that, in general, both strong coupling and anisotropy ef- 
fects play a significant role in the dimensionless ratios, 
and both need to be accounted for. 



V. THE LIMIT OF NEARLY SEPARATE 
BANDS 

When A12 — A21 = 0, there exist two transition tem- 
peratures Tci and Tc2 associated with An and A22, sepa- 
rately, and for several properties, but not all, the super- 
conducting state is the straight sum of the two bands as 
they would be in isolation. Here, we wish to study how 
the integration of the two bands proceeds as A12 and/or 
A21 is switched on. Our first results related to this issue 
are shown in Fig. which has three frames. The top 
frame deals with the normalized specific heat Cs{T)/^T 
as a function of temperature, the middle, the normalized 
inverse square of the penetration depth [Xl{0)/ Xl{T)]'^ 
and the bottom gives the critical field deviation func- 
tion D{t) of Eq. ((23) • In all cases, we have used our 
Lorentzian model for the spectral densities afjF{uj) with 
All = 1 and A22 = 0.5 fixed for all curves. The solid 
curves are for A12 = A21 = 0.0001, short-dashed for 
0.01, and long-dashed for 0.1. In the top two frames. 
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the two separate transitions are easily identified in the 
curves with solid line type. Because of the very small 
value of Ai2 — X21, the composite curve is obviously the 
summation of two subsystems, which are almost com- 
pletely decoupled. However, already for A12 — X21 — 0.01 
which remains very small as compared with the value of 
All and even A22, the second transition (short-dashed 
curve) becomes significantly smeared. The two subsys- 
tems have undergone considerable integration. In par- 
ticular, the second specific heat jump is rounded, be- 
coming more knee-like. Also, the sharp edge or kink in 
the solid curve for the superfluid density is gone in the 
short-dashed curve. Thus, to observe clearly two dis- 
tinct systems, the offdiagonal A's need to be very small. 
Once A12 = A21 =0.1 (long-dashed curve), the integra- 
tion of the two subsystems is very considerable if not 
complete. This does not mean, however, that super- 
conducting properties become identical to those for an 
equivalent one-band system. As long as the afjF{uj) are 
not all the same, there will be anisotropy and this will 
change properties as compared with isotropic Eliashberg 
one-band results. Note that in the solid Eliashberg curve 
of Fig. El a point of inflection remains, as commented on 
by Golubov et al. [iJi] • In the case of the deviation func- 
tion (lower frame), the solid curve shows a sharp cusp 
which is related to the lower transition temperature of 
the decoupled bands but not quite at that value as this 
function is composed from subtracting 1 — (T/Tc)'^ from 
Hc{T)/ Hc{0). However, two distinct pieces of the curve 
exist and notably near Tc the curve has a very different 
curvature from what is normally encountered. In particu- 
lar, the temperature dependence of the solid curve is con- 
cave down at high temperature in contrast to the usual 
case of concave up. As the bands are coupled through 
larger and larger interband A's, the curve moves to a 
shape more consistent with one-band behaviour. How- 
ever, the curve remains negative due to the anisotropy, 
while usually strong coupling would drive the curve pos- 
itive with an overall concave-down curvature^], which is 
illustrated by the dotted curve for which the first band 
dominates, as we describe below. 

The other curves in these figures, dot-dashed and dot- 
ted, are for a = 10 with A12 — 0.1 and A21 — 0.01, and 
a = 0.1 with A12 = 0.01 and A21 = 0.1, respectively. For 
a = 10, the second band with the smaller of the two di- 
agonal values of A has ten times the density of states as 
compared to band 1 with the larger A value. This large 
disparity in density of states can have drastic effects on 
superconducting properties, and further modify both the 
observed temperature dependence and the value of the 
BCS ratios. The second specific heat jump in the dash- 
dotted curve, although smeared, is quite large as com- 
pared with that in the solid or even the dashed curve. 
Also, it is to be noted that beyond the temperature of 
the lower maximum in Cg (T) / jT, the curve shows only 
a very modest increase, refiecting the low value of the 
electronic density of states in band 1, and the ratio of 
the jump at to the normal state is now quite reduced. 



The low density of states in band 1 is also reflected in the 
low value of the penetration depth curve (middle frame, 
dash-dotted curve) in the temperature region above Tc2- 
Finally, we note that while we have chosen a large value 
of a for illustration here, MgB2 has an a = 1.37 which, 
by the above arguments, would tend to accentuate the 
features due to the second band. 
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FIG. 10: Upper frame: Individual contributions from each 
band to the superfluid density [A(0)/Ai,2(r)]^ as a function 
of T/Tco, where Tco is the Tc for the An channel alone, with 
all others zero. Shown are curves for various off'diagonal A21 
with All ~ 1, A22 = 0.5, and A12 — 0.0001. The three pairs 
of curves are for A21 equal to: 0.0001 (solid), 0.01 (short- 
dashed), and 0.1 (long-dashed). The curves which go to zero 
at a lower temperature correspond to [A(0)/A2(r)]^ while 
those which go to zero close to 1 are for [A(0)/Ai(r)]^ Lower 
frame: Now the A21 is held fixed at 0.0001 and the A12 is 
varied. The three pairs of curves are for A12 equal to: 0.0001 
(solid), 0.1 (short-dashed), and 0.2 (long-dashed). Here, the 
ratio of the density of states a has been taken to be 1 for 
convenience of illustrating the curves on the same scale. 

A very different behaviour is obtained when a = 0.1 for 
which case the electronic density of states in the second 
band is reduced by a factor of ten as compared to the flrst 
band. In this case, the dotted curve applies and looks 
much more like a standard one-band case with very sig- 
nificant strong coupling effects AC{Tc)/jTc ~ 2.4. The 
infiuence of band 2 has been greatly reduced. Finally, 
we note that the introduction of the offdiagonal elements 
can change Tc- In particular, the dot-dashed curve ends 
at a considerably reduced value of critical temperature as 
compared with the other curves. This is consistent with 
Fig. ^where we saw that increasing A12 for small values 
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of A21 decreases Tc- On the other hand, for the dotted 
curve for which values of A12 and A21 are interchanged as 
compared to the dash-dotted curve, Tc is hardly affected 
because A12 is small and it is this parameter which affects 
Tc more. The two parameters A12 and A21 do not play 
the same role in Tc or for that matter in the integration 
process of the two bands. This is made clear in Fig. ^| 
which deals only with the penetration depth. What is 
shown are the separate contributions to the superfluid 
density coming from the two bands. In all cases, An = 1 
and A22 = 0.5. In the top frame, A12 = 0.0001 and A21 is 
varied. It is clear that as A21 is increased, the superfluid 
density associated with the second band remains signif- 
icant even above the second transition temperature Tc2 
which is well-defined in the solid curve. This is the oppo- 
site behaviour of what is seen in the lower frame where 
A21 remains at 0.0001 and A12 is increased. In this case, 
Tc changes significantly but the superfluid density associ- 
ated with the second band remains negligible above Tc2. 
Note finally that the relative size of the superfluid den- 
sity in each band will vary with a and vpi, neither of 
which have been properly accounted for in this figure, as 
we wished to illustrate solely the effect of A12 and A21 on 
the issue of integration of the bands and modification of 
Tc 

The changes, with the offdiagonal elements A12 and 
A21, in the temperature dependence of the upper and 
lower gaps are closely correlated with those just described 
for the superfluid density. This is documented in Fig. 1111 
which has two frames. In all cases An = 1 and A22 = 0.5. 
In the top frame, A12 = A21 equal to 0.0001 (sohd), 0.01 
(short-dashed), and 0.1 (long-dashed). The various pairs 
of curves apply to the upper and lower gap ratios. Note 
the long tails in the short-dashed curve (lower gap), still 
small but extending to T = Tc. For the long-dashed 
curve, the lower and upper gaps now have very similar 
temperature dependences, but these are not yet identical 
to BCS. We have already seen in Fig.jTJ for the specific 
case of MgB2, that the lower gap falls below BCS at tem- 
peratures above T/Tc ~ 0.7, which is expected when it 
is viewed as an evolution out of two separate gaps, with 
two Tc values, due to increasing the offdiagonal coupling. 
In the lower frame, we show results for a = 10 (dot- 
dashed) and a = 0.1 (dotted). Again, as expected, the 
two dash-dotted curves show distinct temperature depen- 
dences while for the dotted they are very similar. 

A very similar story emerges when interband impu- 
rity scattering is considered. Results are given in Fig. 1121 
and Fig. El Fig. E| has three frames. Here, An = 1, 
A22 = 0.5, and A12 = A21 = 0.0001, with our Lorentzian 
electron-phonon spectral functions afjF{LLj) described 
previously. The top frame deals with the temperature 
dependence of the normalized superconducting state elec- 
tronic specific heat Cs{T)/jT. The middle frame gives 
the gap ratios of Ai and A2 and thus the curves come in 
pairs, with Ai > A2. And the bottom frame shows the 
deviation function D{t) for the thermodynamic critical 
magnetic field. What is varied in the various curves is 
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FIG. 11: Upper frame: Upper and lower gap ratios, 
2Ai^2/kBTco, versus T/Tco, where Tco is the Tc for the An 
channel alone, with all A's zero. Shown are curves for various 
offdiagonal A's with An = 1 and A22 ~ 0.5. Three pairs of 
curves are for A12 = A21 equal to: 0.0001 (solid), 0.01 (short- 
dashed) and 0.1 (long-dashed). Lower frame: Same as for 
upper frame except now are shown: A12 =0.1 and A21 ~ 0.01 
(i.e., a = 10) (dot-dashed) and A12 — 0.01 and A21 ~ 0.1 
{i.e., a = 0.1) (dotted). 



the interband impurity scattering rate — tti (taken 
to be equal in value, i.e., a = 1). The solid curve, which 
clearly shows two transitions, is for = 0. It is to be 
noted first, that in all cases, offdiagonal impurity scat- 
tering changes the value of the critical temperature, re- 
ducing it to less than 0.8 of its pure value in the case 
of the dot-dashed curve. This decrease in Tc does not 
translate, however, into a steady decrease in the specific 
heat jump at Tc. We see that while the jump initially de- 
creases with increasing — ttn eventually it increases 
and is largest for the dot-dashed curve. Both Watanabe 
and Kita^J and Mishonov et a/.jsj], using only an un- 
renormalized BCS model, find an increase with impurity 
scattering and no initial decrease as is found in the full 
Eliashberg calculation. This is a clear illustration that, 
at minimum, a renormalized BCS formula needs to be 
used to capture the qualitative trend and full Eliashberg 
theory is required if one wishes to be quantitative. It is 
also clear that as interband impurity scattering increases, 
the jump in the specific heat at the second transition, 
seen in the solid curve, is rapidly washed out and little 
remains of this anomaly in the dot-dashed curve. Even 
the long-dashed curve shows little structure in this re- 



16 




-0.06 r 1 

0.0 0.2 0.4 0.6 0.8 1.0 

T/Teo 

FIG. 12: Top frame: Specific heat in the superconduct- 
ing state normalized to the normal state, Cs{T) /'yT , versus 
T/Tco, where Tco is the Tc for the pure case. Here, An = 1, 
A22 = 0.5, and A12 = A21 = 0.0001. Shown are curves for 
varying = ^21 equal to: 0.0 (solid), 0.01 (short-dashed), 
0.2 (long-dashed), and 0.5 (dot-dashed) in units of Tco- No- 
tice that the value of the jump at Tc first dips and then rises 
with impurity scattering. Middle frame: 2/S.\^2 /ksTco versus 
T /Tco- The upper three curves correspond to 2Ai / ksTco and 
the lower three to 2A2/fcsTco, with the curves labelled the 
same way as in the upper frame. Only the first three impu- 
rity cases are shown for clarity. The other progresses in the 
same manner with the Tc reducing further and the gaps mov- 
ing closer to a common value. Bottom frame: The deviation 
function D{t) versus T/Tco, again with the curves labelled the 
same way as in the top frame. 



gion, in analogy to what we found to hold for the case 
of increasing the ofFdiagonal electron-phonon elements. 
Note, however, there remains a point of inflection which 
has moved to higher temperature. Such a shift of the 
inflection point can also be brought about by increasing 
the offdiagonal A's as seen in Fig. |51 

The temperature dependence of the gap ratios (middle 
frame of Fig. I12|) also mirror what we found in Fig. ITTI 
The dashed curves exhibit quite distinct temperature de- 
pendences between Ai and A2 while this is no longer the 
case for the pair of long-dashed curves. Note that, as 
compared to the solid curve, the anisotropy in the gaps 
for the long-dashed curve has been reduced considerably. 
The upper gap has decreased and the lower increased 
even more. The washing out of the gap anisotropy by 
offdiagonal impurity scattering is expected and has been 
studied theoretically |32n53i] and experimentally |5J]. For 



carbon doping, the gaps are seen to merge at about 13% 
for which Tc has been reduced to about 20 K with the 
large gap reducing to its BCS value and the smaller gap 
moving upwards only very little in contrast to our model 
calculations for which the lower gap changes relatively 
more and isotropy is reached at about a 30% reduc- 
tion in Tc- Of course, as one dopes, the electronic den- 
sity of states and the electron-phonon parameters also 
change [ssj. and one needs to include these in addition to 
any interband scattering. 

Finally, the effect of interband ordinary impurity scat- 
tering on the deviation function, shows a behaviour simi- 
lar to that found for paramagnetic impurities in one-band 
superconductors 35J . Initially, as in the other properties, 
the impurities smear the structure related to the second 
transition temperature (in this case the cusp feature in 
the solid curve) and once the two bands are fairly well 
integrated, then like paramagnetic impurities, the effect 
here is to keep the minimum at the same temperature 
but change its value. A key difference though is that, 
in the case of paramagnetic impurities in one-band, the 
extremum in the curve moves from positive (and strong 
coupling) to negative (and weak coupling) because the 
gap is being reduced towards zero. Here, with the two 
bands, the impurities do not reduce T^ and hence the 
gap, to zero, but rather to a finite value related to the 
washing out of the anisotropy between the two bands, 
and hence the extremum in this case will move from neg- 
ative (where it is positioned due to large anisotropy) to 
smaller values, reflecting this. 

Next we turn to the results given in Fig. ^1 which 
shows the temperature dependence of the superfluid den- 
sity for various values of impurity parameters. Again, 
All = 1, A22 = 0.5, and A12 = A21 = 0.0001 with the 
Lorentzian spectra. What is illustrated in these four 
frames is how very different the effect of tj^, tjjj ^12' 
and are. The solid curve is for reference and is the 
pure case. Once again, for the case of varying and 

, we have violated the constraint that their ratio must 
be fixed by the ratio of the density of states. This we can 
do theoretically to decouple and, therefore, illustrate the 
effects of these different scattering channels, but in real 
systems, they would be constrained and the net result 
would be a combination of the effects from both chan- 
nels. The top left frame shows the effect on the superfluid 
density of increasing the impurity scattering in the first 
band (intraband scattering). Such impurities reduce the 
superfluid density in band 1 while leaving band 2 unal- 
tered. In the lower right-hand frame it is the superfluid 
density in the second band that is reduced, leaving the 
first unchanged. Tc is unaffected by intraband impurity 
scattering in isotropic s-wavc superconductors due to An- 
derson's theorem. The top right-hand frame shows that 
increasing reduces the critical temperature as well 
as reduces the superfluid density in band one without, 
however, having much effect on the second band. The 
kink associated with the rise of the second band is hardly 
changed as 1^2 is not the integrating variable, rather it 
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FIG. 13: Effect of impurity scattering on the superfluid den- 
sity [Ao(0)/A(T)]^ plotted versus T/Tco, where T^o and Ao(0) 
are for the pure case. Each frame shows the effect of the 
different type of impurity scattering keeping all other im- 
purity terms equal to zero. The spectrum parameters are: 
All = 1, A22 = 0.5, and A12 = A21 = 0.0001. The solid curve 
in all cases is for the pure case of tf- — Q. In the upper left 
frame: t^^/Tco ~ 0.2 (short-dashed) and 2.0 (long-dashed). In 
the upper right frame: t^j/^^co = 0.2 (short-dashed) and 0.4 
(long-dashed). In the lower left frame: t^^/Tco = 0.02 (short- 
dashed) and 0.1 (long-dashed). In the lower right frame: 
tJj/Tco = 0.2 (short-dashed) and 2.0 (long-dashed). 



is which integrates rapidly the bands as seen in the 
lower left-hand frame. However, in this case, the critical 
temperature is hardly changed and there is little change 
to the curve above T/Tco — 0.7. 
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FIG. 14: Upper frame: Gap ratios for the upper (2Ai/fcsTc) 
and lower gap {2/S.2/kBTc) as a function of A12 for varying A21: 
0.01 (solid), 0.1 (dotted), and 0.3 (dashed). Here An = 0, 
A22 ~ 0. This is for comparison with Suhl et aZ.Q- Lower 
frame: Gap anisotropy, u — A2/A1, versus A12 for the same 
parameters and curve labels as the upper frame. Note that 
A12 > A21 is plotted. With A12 < A21, the roles are simply 
reversed with 1 2 and A2 would become the large gap, etc. 



A22 = 0, Eq. (|12|1 for the coupling A, which determines 
Tc from Eq. Hll|l . simplifies to: 



VI. THE LIMIT OF PURE OFFDIAGONAL 
COUPLING 

While the two-band nature in MgB2, driven by the 
electron-phonon interaction, is well-established, there 
have been many reports of possible two-band supercon- 
ductivity in other systems, including the conventional 
A15 compound NbsSn 56]. With a = 18 K and a 
main gap 2Am ^ 4.9rc, there is specific heat evidence 
for a second gap at O.STc- Other systems are NbSe9[5^. 
Y2C3 and La2C3[53 and possibly a second nonsupercon- 
ducting band in CeCoIn r. 15911 . In the triplet spin state 
superconductor Sr^RuO/tpOll. a small gap is induced in 
the second band. As two-band superconductivity is likely 
to be a widespread phenomenon, not confined to electron- 
phonon systems, it seems appropriate to investigate fur- 
ther an extended range of parameter space for the Xij 's 
and in particular the possibility that the offdiagonal ele- 
ments are the dominant mechanism for superconductiv- 
ity. 

In the limit of pure offdiagonal coupling, where An = 



A: 



(51) 



V A12A21 

and the ratio of the gap to Tc given in Eq. (|21|) becomes 
2Ai 



3.54exp 



A 



u 
A21 



(52) 



The ratio A12/A21 = a* can be taken > 1 and Eq. 
for the gap anisotropy u = A2/A1 written as: 



* 2 1 
a u — 1 — 



A 



dn7 



(53) 



This equation gives u in terms of a* and A. Since by its 
definition 0< u < 1, wlnwis negative so a condition on 
obtaining a solution of H53|l is that 



a'v? — 1 < or a* < 



(54) 



For a trial solution of u — 0.1, this would give 1 < 
a* < 100. For a* = 60, as an example, A = 4.46 and 
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l^xjkBTc — 9.7, which is very large. This occurs for 
Tc/'^\n ^ 10~^, using ln(1.13a;in/Tc) = A, which is in 
the weak coupHng regime. However, to achieve an up- 
per gap ratio value greater than 11 or so, will correspond 
an unrealistically small value of Tc/uJin (of order lO^^*', 
for example). In Fig. 1141 we show results in the upper 
frame for 2Ai_2/fcsT'c versus A12 for various A21 values. 
In the lower frame, we show u versus A12. The differ- 
ence between Fig. ITll and Fig. 13 shows that large values 
of 2Ai/fcBT'c are more naturally obtained in the pure 
offdiagonal regime and are associated as well with small 
values of u and the weak coupling regime. This latter fea- 
ture implies that there will be no further strong coupling 
corrections to an already large gap ratio. We have also 
calculated the thermodynamics and superfluid density in 
this regime, for a range of parameters, but have found 
these properties to show quite ordinary behaviour and 
have discovered no new physics. For the sake of brevity, 
we present none of these results but instead note that in 
this limit the x's are 



X2 



A21 
A12 
A12 
A21 



1 

a* 



= 2(1 + "*)' 



(55) 
(56) 



with a* — A12/A21, and hence, the various dimensionless 
ratios are: 



AC 

Wc 



= 1.43 



4a* 



(1 + a*)^ 



and 



and 



hciO) = 



Tr TT 



All /7C(3) 
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(1 + a*) 



v*2 



\y'LiTc)\Tc 



1 (l + a/?)(l + a*) 
4 a* + Pa 



(57) 



(58) 



(59) 



where /3 = Vp2i^ + ^i2)/vpi{l + A21). The ratio for the 
zero temperature critical field of Eq. H38|l does not change 
its form and so is not repeated here. These ratios behave, 
qualitatively, no differently from what we found in section 
III. A difference worth noting is the following. In linear 
order, the effect of interband impurity scattering on T^, 
takes the form (|16ll8f) : 



AT, 

Tco 




(60) 



which is always negative and larger for paramagnetic 
than for normal impurities. It can also be very large 
for 3> 1. This is another distinction between pure 
offdiagonal coupling and MgB2, for example. In obtain- 
ing H60|l . we have used the fact that A12/A21 = P12/P21 
and A12/A21 > 1. 



VII. CONCLUSIONS 

We have calculated thermodynamics, gap anisotropy 
and penetration depth for a two-band Eliashberg super- 
conductor. For the parameters appropriate to MgB2 
which are obtained from first principle band structure 
calculations of the electron-phonon spectral functions, 
we find good agreement with the existing experimental 
data. We reduce the Eliashberg equations to a renor- 
malized BCS form by application of the two-square-well 
approximation. Comparison of these results with those 
from the full Eliashberg equations allows us to determine 
strong coupling corrections, which we find to be signif- 
icant in MgB2. When the parameters for the electron- 
phonon interaction are moved away from those specific to 
MgB2, the strong coupling corrections can become much 
larger, and superconducting properties reflect this fact, 
as well as the change in anisotropy between the bands. 
Within the A^* approximation, we derive simple ana- 
lytic expressions for the various dimensionless BCS ratios 
which would be universal in the one-band case, but are 
not in the two-band one. They depend on the anisotropy 
and particularly on the ratio of the electronic density of 
states in the two bands. The anisotropy in the ratio of 
the two gaps at zero temperature is investigated and is 
found to increase as A22 is reduced and made repulsive, in 
which case the existence of superconductivity in the first 
band, and the offdiagonal coupling to it, induces a gap in 
a band which would, on its own, not be superconducting. 

We have paid particular attention to the limit of nearly 
decoupled bands, i.e., small interband coupling, with the 
superconductivity originating from An and A22 in the 
first and second band, respectively. When Ai2,A2i 0, 
there are two transitions at Td and Tc2 and two specific 
heat jumps. As the interband coupling is turned on, the 
two bands become integrated and the second transition 
smears. We have found that the two parameters, A12 and 
A21 , have very different effects on the smearing of the sec- 
ond transition and on T,. A12 largely modifies T,, reduc- 
ing it, whereas, A21 alters the lower temperature region 
around the second transition. Only very small values of 
A21, as compared with An and A22, are needed to cause 
large changes in the region around Tc2. It was found that 
a small amount of interband impurity scattering can also 
significantly smear the second transition, and so reduce 
the distinction between the two bands. However, even 
when the two bands are well-integrated and a sharp sec- 
ond transition is no longer easily discernible, this does 
not imply that the superconducting properties become 
those of a one-band superconductor. Anisotropy remains 
and this affects properties. 

In view of the possible widespread occurrence of two- 
band superconductivity, even for systems with exotic 
mechanism not necessarily due to the electron-phonon in- 
teraction, we deemed it of interest to consider the case of 
zero intraband coupling. An — A22 = 0, with supercon- 
ductivity due only to the interband A12 and A21, which 
need not have the same value. When these are very differ- 
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ent, the resulting gaps are quite different from each other 
and the ratio of Ai to can become large particularly in 
the weak coupling limit. This is a distinguishing feature 
of pure offdiagonal coupling. Another distinguishing fea- 
ture is the possibility of a rapid reduction of Tc towards 
zero by interband impurity scattering, as compared with 
the case for which the diagonal elements play the leading 
role. 
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